Pseudo-magnetic catalysis of the time-reversal-symmetry breaking in graphene 
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A finite flux of the (time-reversal-symmetric) pseudo-magnetic field, which would represent the 
effect of a bulge in the graphene sheet for example, is shown to be a catalyst for spontaneous 
breaking of the time-reversal symmetry of Dirac fermions in two dimensions. Possible experimental 
consequences of this effect for graphene are discussed. 



I. INTRODUCTION 

It is well appreciated that the graphene sheet provides 
a particularly simple and physically relevant table-top 
realization of the two-dimensional (pseudo) relativistic 
electron dynamics I']. In particular, the Dirac nature of 
graphene's quasiparticles provides these low-energy exci- 
tations with an extra protection from the usual effects of 
electron-electron interactions. The semi-metallic, non- 
interacting ground state of electrons in graphene may 
be understood as a Gaussian fixed point in the space 
of coupling constants, stable in all directions. Never- 
theless, the ground state can in principle be turned into 
one with a broken symmetry at a finite, and, relative to 
the bandwidth, typically large interaction Q]. This way, 
for example, the system would acquire a finite staggered 
density, or a staggered magnetization, at a large nearest- 
neighbor and on-site repulsion, respectively. Both of 
these order parameters correspond to finite "masses" of 
the Dirac fermions that reduce the chiral ("valley", or 
"pseudo-spin") SU{2) symmetry of the linearized Hamil- 
tonian down to U{1). The two-dimensional Dirac Hamil- 
tonian, however, admits an additional mass-term that 
is invariant under the chiral symmetry, but odd under 
time-reversal Q. It has been argued recently that such 
a time-reversal-symmetry-breaking mass would be gen- 
erated dynamically at a large second-nearest-neighbor 
repulsion between electrons on a honeycomb lattice [1]. 
The type of mass, or an order parameter, that would 
eventually open up at strong coupling, seems to depend 
therefore on the non-universal details of the interactions 
on the atomic scale. 

Increasing the density of low-energy states is expected 
to enhance the effects of interactions on Dirac fermions. 
A manifestation of this general principle is the "mag- 
netic catalysis" , by which the chiral-symmetry-breaking 
mass is induced at an infinitesimal favorable interaction 
in a uniform magnetic field Q i ■ This mechanism is at 
the heart of several recent theories of some of the quan- 
tum Hall eflFects observed in graphene 0, S Q- Mag- 
netic field cannot catalyze the time-reversal-symmetry- 
breaking mass, however [loj. The purpose of this work 
is to show that the fiux of the (non-abelian) pseudo- 
magnetic field plays the role of such a catalyzer. I demon- 
strate that in the presence of a finite flux of the non- 
abelian gauge field an infinitesimal favorable interaction 
would lead to the spontaneous breaking of the time- 
reversal symmetry of the ground state of two-dimensional 



Dirac fermions. This result is a general property of Dirac 
fermions in two dimensions, and as such it is indepen- 
dent of the specific nature of the underlying physical 
system. Nevertheless, its immediate significance derives 
from the notion that a component of such a pseudo- 
magnetic field represents the effect of smooth height vari- 
ations of graphene's surface on the electron dynamics 
With this possible application in mind I consider both 
the limits of a perfectly uniform and a spatially local- 
ized pseudo-magnetic flux. It is found that even the lat- 
ter catalyzes a finite, but local, time-reversal-symmetry- 
breaking mass. Experimental conditions for this non- 
intuitive manifestation of the coupling between the elec- 
tronic and mechanical degrees of freedom in graphene are 
discussed. 



II. DIRAC HAMILTONIAN AND THE 
TIME-REVERSAL OPERATOR FOR GRAPHENE 

Let me establish the notation first. Consider the Dirac 
Hamiltonian for the four-component massless fermions in 
two spatial dimensions: 

H[A'>,A]=^Jol^iP^'A'^-A,), (1) 

where the repeated index i = 1,2 is summed over, and 
Ai is the general non-abelian SU{2) gauge field 

(2) 



where 735 = 17375. A^ is the C/(l) (abelian) compo- 
nent that represents the physical magnetic field, whereas 
Aj, j = 3,5,35 multiply the three generators of the 
chiral SU (2) symmetry 11] of the free Dirac Hamilto- 
nian H[0, 0]. The five gamma-matrices satisfy {7^, j^} = 
2(5^1/, /i = 0, 1, 2, 3, 5, and we will define them here to be 
all Hermitian. In our units, h, = e = c = \. 

The general mass-term that can be added to the Hamil- 
tonian (1) which violates the SU{2) chiral symmetry is 

given by M = m ■ M, where AI = (70, «7o73, *7o75) is a 
vector under the chiral transformations. An additional 
mass-term may then be defined to be a chiral scalar: 
toM, with M = 17172. It is easy to check that the set of 
all linearly independent matrices that anticommute with 
the free Dirac Hamiltonian H[0, 0] is exhausted by M and 
M, which therefore represent all the possible mass-terms. 

An important role in the discussion will be played by 
the time-reversal symmetry of the free Dirac Hamilto- 
nian. As usual, the time-reversal is represented by an 
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anti-unitary operator /( = UtK, where Ut is unitary, and 
K stands for the complex conjugation [l^. Although 
everything that will be discussed hereafter will be man- 
ifestly representation- independent, to exhibit the time- 
reversal operator one needs some representation of the 
7-matrices. We prefer the "graphene representation" in- 
troduced earlier 0], in which 70 = l2®o'z, 71 = (Jz® (Jy, 
= h ® cTx, 73 = tJa; (8) ay, and 75 = cr^ (g) ay, with 
{hi a} as the standard Pauli basis in the space of two- 
dimensional matrices. In this representation the time- 
reversal invariance of the free Dirac Hamiltonian i7[0, 0] 
and of the general chiral-symmetry-breaking mass M de- 
termines the unitary part of the time-reversal operator 
uniquely to be: 

Ut = i7i75 = {ax ® h)- (3) 

Postulating time-reversal invariance of both i/[0, 0] and 
M is motivated by the fact that these operators represent 
the low-energy limit of a completely real lattice Hamil- 
tonian [13. As an immediate consequence, the chiral- 
symmetry preserving mass M must be odd under time- 
reversal. This then is also in accord with the concrete 
lattice realization Q of M. 

III. MAGNETIC CATALYSIS OF CHIRAL 
SYMMETRY BREAKING 

We begin by reformulating the mechanism of the mag- 
netic catalysis in purely algebraic terms. Consider the 
Hamiftonian iJ[A°,0] with A° ^ 0. By virtue of repre- 
senting the physical magnetic field H[A^ , 0] has the time- 
reversal symmetry broken, but the chiral symmetry pre- 
served. In general, the spectrum of H[A^ , 0] will contain 
states with exactly zero energy T^]. Let us denote that 
zero-energy subspace of the full Hilbert space Ho- Ho 
is invariant under the generators of the chiral symme- 
try which by definition all commute with ff[j4°,0], but 
also under the operators that anticommute with H[A^ , 0] , 
such as M and M. If we denote the trace of an operator 
within Tio as Ttq, it follows that 

TroM = 0. (4) 

This is because for each component of M there exists 
an operator which leaves Tio invariant and anticommutes 
with it In the basis of Ho which diagonalizes a 

chosen component of Af the number of states with the 
eigenvalue -1-1 is thus equal to the number of those with 
the eigenvalue —1. Since one can write the ground state 
expectation value of a traceless operator that anticom- 
mutes with the Hamiltonian, such as M, as [l^, [l3|, 

^^) = ^[E - E ]<„(^)M<fo,„(x), (5) 

n.occup n, empty 

with {"I>o,n} a-s a basis in Hq, we see that occupying all 
the zero-energy eigenstates and leaving the —1 eigen- 
states empty creates the maximal spatial average of the 



above order parameter. At half-filling and in the non- 
interacting system, of course, the ground state is highly 
degenerate, and averaging over all the ground states ul- 
timately leads to vanishing order. Nevertheless, in pres- 
ence of even an infinitesimal interaction that favors a 
finite (A/), the non- interacting ground state is unstable 
towards a new non-degenerate ground state with all -f 1 
states shifted slightly downward and all —1 states up- 
ward in energy, so that the chiral symmetry would be- 
come spontaneously broken. 

In a uniform magnetic field the above mechanism leads 
to a constant chiral-symmetry-breaking order parameter, 
and a gap in the spectrum at an infinitesimal favorable in- 
teraction between Dirac fermions, i. e. "magnetic catal- 
ysis" ^ , [T^ . Obviously the mechanism is quite general, 
and as will be discussed here it will be operative even if 
the magnetic field is not uniform, as long there is a finite 
support of the energy spectrum at zero. 

Before turning to our main subject, it is instructive 
to see why the above mechanism does not lead to the 
catalysis of the chirally symmetric order parameter (M) . 
First, note that unlike M, M commutes with allihe other 
operators that leave Ho invariant, i. e. the generators of 
SU{2) and M, so it does not readily follow that its trace 
within Ho must vanish. In fact, since 

H^[A", 0] = (ft - A^f + Me,,d,Al (6) 

at least for an uniform (and say, positive) magnetic field 
it is obvious that all states in Ho have the same (— 1) 
eigenvalue of AI. That this is generally true may be seen 
by rewriting the Dirac Hamiltonian in the magnetic field 
and in the Coulomb gauge diA^ — as 

H[A°, 0] = e-'^(^)^^i7[0, 0]e-^(^)*'^, (7) 

where = eijdjX- This (non- unitary) transformation 
tells us that the zero-energy states of H[A''\ 0] and of the 
free Hamiltonian are related as 

a>o,„[^°](f)oce^(-^^^$o.„[0](f). (8) 

Since for a total flux F (in units of hc/e) localized near 
the origin, at large \x\ x{^) — ^^Inlafl, the last equation 
implies that only the zero-energy eigenstates of i/[0, 0] 
with the eigenvalue —1 of M may lead to normalizable 
states of i?[v4°,0]. All the states in Ho arc thus the 
— 1 eigenstates of M, even for an arbitrary configura- 
tion of the magnetic field. Eq. (5) then implies that 
J dx{M) = at half- filling for any occupation of the 
zero-energy states. 

To summarize, at the filling one-half, the ground state 
of the Dirac Hamiltonian i/[A°,0] in Eq. (1) in pres- 
ence of a finite magnetic flux, which breaks the time- 
reversal and preserves the chiral symmetry, is inherently 
unstable towards the dynamical generation of the mass, 
that would break the chiral while preserving the time- 
reversal symmetry. I show next that when the physi- 
cal (abelian) magnetic field vanishes and only the (non- 
abelian) pseudo- magnetic field is present, the same may 
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be said, only with the "time-reversal" and the "chiral 
symmetry" in the last sentence exchanged. 



IV. PSEUDO-MAGNETIC CATALYSIS OF 
TIME-REVERSAL SYMMETRY BREAKING 

The time-reversal symmetry, being broken by the mag- 
netic field, did not play any role in the above discussion 
of the usual magnetic catalysis. Let us consider now the 
Hamiltonian H[0,A] in Eq. (1), with A ^ 0. Since the 
time-reversal operator /( anticommutes with all the gen- 
erators of SU{2), it follows that i/[0,A] is even under 
time-reversal. For a general non-abelian gauge configu- 
ration A the chiral SU{2) symmetry of the free Hamil- 
tonian will be completely broken, and if A is everywhere 
proportional to one and the same linear combination of 
the generators, it will be reduced to U{1). M, however, 
still always anticommutes with H[0,A]. Ti.Q in this case 
will thus still be invariant under M, as well as under 
It- As these two operators anticommute, it immediately 
follows that when ^" = and A, ^ 0, 



TroM = 0. 



(9) 



Substituting M for M in Eq. (5), it is now the chiral- 
symmetry-preserving, time-reversal-symmetry-breaking 
order parameter (M) that will become catalyzed in pres- 
ence of an infinitesimal favorable interaction. 

Chiral-symmetry-breaking mass, in turn, is not cat- 
alyzed by the non-abelian gauge field. Assume for exam- 
ple that Ai = Ap735, and A° = 0. Similar to Eq. (7) 
one can then write [l^ 



(10) 



where now Af^ — Eijdjcj). In analogy with the Eq. (8) 
it follows that all the states in Ho now have the same 
eigenvalue of 70 , and it is the chiral-symmetry-breaking 
order parameter J dxljo) that vanishes at half-filling. 

To see the dynamical consequences of the above alge- 
bra more explicitly, consider the Lagrangian density 

C = ^{x)^^{d^~Afj,,)^ix)~-^i¥{x)M^ix)f (11) 

with an interaction g > 0, fi — 0,1,2, x = {xo,x), xq 
as the imaginary time, and the quenched component of 
the non-abelian gauge field Af^{x) ^ 0. Introducing the 
Hubbard-Stratonovich field this can be rewritten as 



C^^ix)^^id^-Af{x)js5)^ix) + 

-^m^(x) - rh{x)'^\x)M'i/{x). 
^9 



(12) 



The mean-field theory of the above interacting problem 
would amount to minimization of the corresponding ac- 
tion J Cdx with respect to to(x), or equivalently, to de- 
termining the ground-state expectation value 



(^'t(a;)M*(2;)) = {m{x))/g, 



(13) 



self-consistently. An uniform {m{x)) may be understood 
as the time-reversal symmetry breaking order parameter 
of ref. [3]. For a constant pseudo-magnetic field B"^^ — 
d\AS^ —diA^ Avi full analogy with the standard magnetic 
catalysis 0], [l] we then find 



(14) 



where the first term derives from the split zero-energy 
level, and the term 0{g) is due the other Landau levels. 
For an inhomogeneous B^^ (x) the self-consistent calcula- 
tion can be performed only numerically. Here we circum- 
vent this hurdle by dropping the self-consistency require- 
ment and minimizing the action with respect to only a 
uniform in. This may be understood as a variational 
calculation, or as the exact solution of the Berlin-Kac 
version of the theory [l^, in which the contact interac- 
tion in Eq. (11) is replaced with the interaction of an 
infinite range |20|: 



_9_ 

2Vl 



dy{^\x)M^{x)){^\y)M^{y)), (15) 



with as the area of the system. The uniform ansatz 
becomes an exact solution of the modified theory in the 
thermodynamic limit 17 cxd. In either case there is a 
gap of 2m in the spectrum, which satisfies 



m 
9 



F 



(g2+ ^,2)1/2' 



(16) 



with JV{e) as the exact density of states per unit area 
of the non-interacting Dirac fermions in the fiux of Af^, 
at e ^ 0. Since at low energies we expect that A/'(e) oc 
£(2-^)/^ 1^21], for z < 2 the second term may be neglected 
at a weak coupling, and m is finite in the thermodynamic 
limit only in the case of an extensive flux, F (x fl 
Nevertheless, even if F is finite the expectation value of 
the time-reversal symmetry breaking order parameter is 
finite and equal to 

hm {^\x)M^ix)) = IJ2 <„(^)*o,„(f). (17) 

■Ho 

A finite pseudo-flux selects the time-reversal-symmetry 
broken ground state out of the degenerate manifold, in 
close parallel with the standard formalism of spontaneous 
symmetry breaking in statistical physics [23| . 

To illustrate the local character of the order parameter 
for finite F, let us exhibit the sum in the last equation 
for the particular pseudo-magnetic field 



B-'^r) 



2F 



i?2(l + (r/i?)2)2 



(18) 



In the graphene representation the zero-energy state with 
±1 eigenvalue of M are then $1 _(x) = f*{x){0, 1, 0, 0), 
and $|'j^(f) = /„(£)(0, 0,0,1), where 



7r-li?-2("+l)(2; + iy)» 



^(3{n + l,F-n- 1)(1 + {r/Rf) 



2\F ' 



(19) 
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with the integer n < F. Note that the $n,+ = h^n.-- 
For an integer flux F then the sum in Eq. (17) can be 
exactly performed with the result 

lim {¥{x)M^{x)) = (1 - F-^)B^^{r). (20) 

O — >oo 

For a general localized flux the precise proportionality 
between the order parameter and the field obtains only 
in the limit F ^ 1 24]. The order parameter, however, 
is always localized in the region of fiux. 

V. EXPERIMENTAL CONSEQUENCES 

Finally, let us address possible consequences of the 
above results for graphene. As mentioned in the in- 
troduction, the time-reversal-symmetry-breaking mass 
is favored by the second- nearest-neighbor repulsion 
whereas the competing chiral-symmetry-breaking masses 
are preferred by the nearest-neighbor repulsion between 
electrons. With the electron spin included, chiral- 
symmetry-breaking mass with the opposite sign for the 
two spin components, which corresponds to staggered 
magnetization, is also preferred by the, most likely the 
strongest, on-site repulsion As one has little control 
over the size of the interaction couplings and can hope 
only to alter the bandwidth, the possible instability to- 
wards the time-reversal-symmetry-breaking mass with- 
out any gauge fields seems likely to be inferior to the one 
towards chiral-symmetry breaking. An "application" of 
the pseudo-magnetic fiux, however, changes this, since it 
is only the time-reversal-symmetry-breaking mass that is 
catalyzed by it at weak interactions. 

A crude estimate of the locally catalyzed gap gives 
m « VB^^ / Biatt, where Biatt ~ lO'^T is the characteris- 
tic lattice magnetic field scale, and V is the strength of 
the second-nearest-neighbor repulsion. A (single) wrin- 



kle which tends to spontaneously form on agraphene 
sheet would already lead to B^^ ~ IT [11], [20], so to- 
gether with an estimate of ~ (1 — 5)eV [3, m ~ 
(0.1 — 0.5)'meV. A randomly wrinkled graphene corre- 
sponds to zero total flux, of course, and so {rh{x)) — 0. 
To produce a flnite net pseudo-magnetic fiux one needs 
to deliberately bulge the graphene sheet, which accord- 
ing to the above estimate should push the gap well into 
the meV range. 

The pseudo-magnetic catalysis described here is stable 
with respect to deviations from half-filling, i. e. for the 
chemical potential smaller than the generated mass. 



VI. SUMMARY 

To conclude, I described the mechanism complemen- 
tary to the usual magnetic catalysis: a finite net fiux of 
a component of the non-abelian gauge field, which pre- 
serves the time-reversal and breaks the chiral symmetry 
of the free Dirac Hamiltonian, serves as a catalyst of 
the time-reversal-symmetry-breaking, chiral-symmetry- 
preserving order parameter. This could lead to local 
spontaneous breaking of the time reversal symmetry in 
graphene where such a pseudo-magnetic field is provided 
by a bulge in graphene's plane, due to the second-nearest- 
neighbor repulsion term in the lattice Hamiltonian. The 
magnitude of the effect should be large enough for the 
gap in the local density of states to become observable 
by scanning tunneling microscopy, for example [27| . 
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